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A huge amount of good quality astrophysical data converges towards the picture of a spatially 
flat universe undergoing the today observed phase of accelerated expansion. This new observational 
trend is commonly addressed as Precision Cosmology. Despite of the quality of astrophysical surveys, 
the nature of dark energy dominating the matter-energy content of the universe is still unknown and 
a lot of different scenarios are viable candidates to explain cosmic acceleration. Methods to test these 
' cosmological models are based on distance measurements and lookback time toward astronomical 

, objects used as standard candles. 1 discuss the characterizing parameters and constraints of three 

, different classes of dark energy models pointing out the related degeneracy problem which is the 

signal that more data at low (z ~ 1), medium (1 < z < 10) and high (10 < z < 1000) redshift 
are needed to definitively select realistic models. 
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I. INTRODUCTION 



The increasing bulk of data that have been accumulated in the last few years have paved the way to the emergence 
of a new standard cosmological model usually referred to as the concordance model. The Hubble diagram of Type 
la Supernovae (hereafter SNela), measured by both the Supernova Cosmology Project>i and the High-z Team^ up 
O ' to redshift z ~ 1, was the first evidence that the universe is undergoing a phase of accelerated expansion. On the 
^ . other hand, balloon born experiments, such as BOOMERanG* and MAXIMA^, determined the location of the first 
^ ■ and second peak in the anisotropy spectrum of cosmic microwave background radiation (CMBR) strongly pointing 
out that the geometry of the universe is spatially flat. If combined with constraints coming from galaxy clusters on 
the matter density parameter ^Im, these data indicate that the universe is dominated by a non-clustered fluid with 
negative pressure, generically dubbed dark energy, which is able to drive the accelerated expansion. This picture has 
^ , been further strengthened by the more precise measurements of the CMBR spectrum, due to the WMAP experiment^, 
and by the extension of the SNela Hubble diagram to redshifts higher than 1~. After these observational evidences, an 
overwhelming flood of papers, presenting a great variety of models trying to explain this phenomenon, has appeared; 
' in any case, the simplest explanation is claiming for the well known cosmological constant A'*. Although the best fit 
to most of the available astrophysical datcfc^, the ACDM model failed in explaining why the inferred value of A is so 
tiny (120 orders of magnitude lower) compared to the typical vacuum energy values predicted by particle physics and 
J — ' why its energy density is today comparable to the matter density (the so called coincidence problem). As a tentative 
, solution, many authors have replaced the cosmological constant with a scalar field rolling down its potential and giving 
rise to the model now referred to as quintessence. Even if successful in fitting the data, the quintessence approach to 
dark energy is still plagued by the coincidence problem since the dark energy and matter densities evolve differently and 
reach comparable values for a very limited portion of the universe evolution coinciding at present era. In this case, the 
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coincidence problem is replaced with a fine-tuning problem. Moreover, it is not clear where this scalar field originates 
. . . from thus leaving a great uncertainty on the choice of the scalar field potential. The subtle and elusive nature of dark 
energy has led many authors to look for completely different scenarios able to give a quintessential behavior without 
the need of exotic components. To this aim, it is worth stressing that the acceleration of the universe only claims 
for a negative pressure dominant component, but does not tell anything about the nature and the number of cosmic 
fiuids filling the universe. This consideration suggests that it could be possible to explain the accelerated expansion 
by introducing a single cosmic fluid with an equation of state causing it to act like dark matter at high densities 
and dark energy at low densities. An attractive feature of these models, usually referred to as Unified Dark Energy 
(UDE) or Unified Dark Matter (UDM) models, is that such an approach naturally solves, al least phenomenologically, 
the coincidence problem. Some interesting examples are the generalized Chaplygin gasi^, the tachyon fieldii and 
the condensate cosmologyiS,. A different class of UDE models has been proposed^'^ where a single fiuid is considered 
whose energy density scales with the redshift in such a way that the radiation dominated era, the matter domination 
era and the accelerating phase can be naturally achieved. It is worth noting that such class of models are extremely 
versatile since they can be interpreted both in the framework of UDE models and as a two-fluid scenario with dark 
matter and scalar field dark energy. The main ingredient of the approach is that a generalized equation of state can 
be always obtained and observational data can be fitted. Actually, there is still a different way to face the problem of 
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cosmic acceleration. As stressed in Lue et al^, it is possible that the observed acceleration is not the manifestation 
of another ingredient in the cosmic pie, but rather the first signal of a breakdown of our understanding of the laws 
of gravitation. From this point of view, it is thus tempting to modify the Friedmann equations to see whether it is 
possible to fit the astrophysical data with models comprising only the standard matter. Interesting examples of this 
kind are the Cardassian expansion^^ and the DGP gravity^-. Moving in this same framework, it is possible to find 
alternative schemes where a quintessential behavior is obtained by taking into account effective models coming from 
fundamental physics giving rise to generalized or higher order gravity actions^''. For instance, a cosmological constant 
term may be recovered as a consequence of a non- vanishing torsion field thus leading to a model which is consistent 
with both SNela Hubble diagram and Sunyaev- Zel'dovich data coming from clusters of galaxies —. SNela data could 
also be efficiently fitted including higher-order curvature invariants in the gravity Lagrangian^-. It is worth noting that 
these alternative schemes provide naturally a cosmological component with negative pressure whose origin is related 
to the geometry of the universe thus overcoming the problems linked to the physical significance of the scalar field. It 
is evident, from this short overview, the high number of cosmological models which are viable candidates to explain 
the observed accelerated expansion. This abundance of models is from one hand the signal of the fact that we have 
a limited number of cosmological tests to discriminate among rival theories, and from the other hand that a urgent 
degeneracy problem has to be faced. To this aim, it is useful to remember that both the SNela Hubble diagram and 
the angular size - redshift relation of compact radio sources^ are distance based methods to probe cosmological models 
so then systematic errors and biases could be iterated. From this point of view, it is interesting to look for tests based 
on time-dependent observables. For example, one can take into account the lookback time to distant objects since 
this quantity can discriminate among different cosmological models. The lookback time is observationally estimated 
as the difference between the present day age of the universe and the age of a given object at redshift z. Such an 
estimate is possible if the object is a galaxy observed in more than one photometric band since its color is determined 
by its age as a consequence of stellar evolution. It is thus possible to get an estimate of the galaxy age by measuring 
its magnitude in different bands and then using stellar evolutionary codes to choose the model that reproduces the 
observed colors at best. A quite similar approach was pursued by Lima & Alcania^i who used the age (rather than 
the lookback time) of old high redshift galaxies to constrain the dark energy equation of state (see alsc^^). The same 
method was then applied also to braneworld models^"^ and the Chaplygin gas^^. It is worth noting, however, that the 
estimate of the age of a single galaxy may be affected by systematic errors which are difficult to control. Actually, 
this problem can be overcome by considering a sample of galaxies belonging to the same cluster. In this way, by 
averaging the estimates of all galaxies, one obtains an estimate of the cluster age and reduces the systematic errors. 
Such a method was first proposed by Dalai et ali^ and then used by Ferreras et ali^ to test a class of models where a 
scalar field is coupled with the matter term giving rise to a particular quintessence scheme. These lectures are devoted 
to review different classes of dark energy models discussing some methods to constrain them toward observational 
data. Far from being exhaustive and complete, my aim is to point out the degeneracy problem and the fact that we 
need further and self-consistent observational surveys at all redshifts to remove it. The layout of these notes are the 
following. In Sect. 2, I briefly present three classes of cosmological models pointing out the main quantities which 
we need for distance measurements and lookback time tests. Methods to fit and constrain models are discussed in 
Sect. 3. In Sect. 4, I take into account a sample of galaxy clusters as an example of data useful to fit dark energy 
models at low redshifts. As exercise, in Sect. 5, the above classes of models are fitted against the sample of galaxy 
clusters considering also SNela and WMAP data and results. The aim of this section is pedagogical: the goal is to 
show the emergence of the degeneracy problem. Summary and conclusions are given in Sect. 6. 



II. DARK ENERGY MODELS 



As discussed in the Introduction, many rival theories have been advocated to fit the observed accelerated expansion 
and to solve the puzzle of the nature of the dark energy. As a simple classification scheme , we may divide the different 
cosmological models in three wide classes. According to the models of the first class, the dark energy is a new ingredient 
of the cosmic Hubble flow, the simplest case being the ACDM scenario and its quintessential generalization which we 
will refer to as QCDM models. This is in sharp contrast with the assumption of UDE models (the second class) where 
there is a single fluid described by an equation of state comprehensive of all regimes of cosmic evolution-'^^ which I 
will consider here referring to it as the parametric density models or generalized EoS^'^ models. Finally, according to 
the third class models, accelerated expansion is the first evidence of a breakdown of the Einstein General Relativity 
(and thus of the Friedmann equations) which has to be considered as a particular case of a more general theory of 
gravity. As an example of this kind of models, we will consider the f{R) - gravit y^^'^^ . Far from being exhaustive, 
considering these three classes of models allow to explore very different scenarios proposed to explain the observed 
cosmic acceleration. In the following we will sketch these three dark energy approaches and derive some of the main 
quantities which we need for matching the observations. 
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FIG. 1: The CMBR anisotropy spectrum for different values of w. Data points are the WMAP measurements and the best fit 
is obtained for ly ~ — 1. 

A. The ACDM model and its generalization toward quintessence (QCDM models) 

Cosmological constant A has become a textbook candidate to drive the accelerated expansion of the spatially flat 
universe. Despite its conceptual problems, the ACDM model turns out to be the best fit to a combined analysis of 
completely different astrophysical data ranging from SNela to CMBR anisotropy spectrum and galaxy clustering^i^li^. 
As a simple generalization, one may consider the QCDM scenario in which the barotropic factor w = p/ p takes at a 
certain epoch a negative value with w — —1 corresponding to the standard cosmological constant. Testing whether 
such a barotropic factor deviate or not from —1 is one of the main issue of modern observational cosmology. How 
such a negative pressure fluid drives the cosmic acceleration may be easily understood by looking at the Friedmann 
equations : 



{PM + Pq) , 



(1) 



(2) 



where the dot denotes differentiation with respect to cosmic time H is the Hubble parameter and the universe 
is assumed spatially flat as suggested by the position of the first peak in the CMBR anisotropy spectrum (see also 
Figd]). 

From the continuity equation, p + 3H{p + p) =0, we get for the i-th fluid with pi = WiPi : 



(3) 



where z=l/a— lis the redshift, Qi = Pi/pcrit is the density parameter for the i-th fluid in terms of the critical 
density Pcrit — 01^/8^0 and, hereafter, I label all the quantities evaluated today with a subscript 0. Inserting this 
result into Eq.([T]), one gets : 



H{z) = Ho^^Mflil + zY + r!Q,o(l + z)3(i+») 



(4) 
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Using Eqs.(IT]), ([2]) and the definition of the deceleration parameter q = —aa/a?, one finds : 



1 3 



(5) 



The SNela Hubble diagram, the large scale galaxy clustering and the CMBR anisotropy spectrum can all be fitted 
by the ACDM model with {yiM,oi^Q) — (0.3,0.7) thus giving go — —0.55, i.e. the universe turns out to be in an 
accelerated expansion phase. The simplicity of the model and its capability of fitting the most of the data are the 
reasons why the ACDM scenario is the leading candidate to explain the dark energy cosmology. Nonetheless, its 
generalization, QCDM models, i.e. mechanisms allowing the evolution of A from the past are invoked to remove the 
A-problem and the coincidence problem. 

B. Generalizing the EoS: the parametric density models 

In the framework of UDE models, it has been proposed^"^ a phenomcnological class of models by introducing a 
single fluid^ with energy density : 



p{a) = Art 
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(6) 



with < a < /3, s and b (with s < b) two scaling factors and Anorm a normalization constant. For several applications, 
it is useful to rewrite the energy density as a function of the redshift z. Replacing a — {1 + z)~^ in Eq.®, we get : 



p{z) = A„ 
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having defined = 1/s — 1 and = 1/6 — 1. It is easy to show that p cx a^'^ for a << s, p oc a~" for s << a << b 
and p cx const, for a » b. By setting (a,/?) — (3,4) the energy density smoothly interpolates from a radiation 
dominated phase to a matter dominated period finally approaching a de Sitter phase. The normalization constant 
Anorm may be estimated by inserting Eq. ^ into Eq. ^ and evaluating the result today. This gives : 
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Pcritfi 



(1 + s)'3-" (1 + fe") 



(8) 



The continuity equation may be recast in a form that allows to compute the pressure and the barotropic factor 
w = pI p obtaining— : 



w 



[{a - 3)a +{13 ~ 3)s]6" - [3(a + s) + (a - /3)s]a" 
3(a + s)(a" +6") 



(9) 



which shows that the barotropic factor strongly depends on the scale factor (and hence on the redshift). Combining 
the Friedmann equations, the deceleration parameter q = (1 + 3u')/2 is 



_ [{a - 2)a + (/3 - 2)s]b°' - [2{a + s) + (a - /3)s]a" ^ 
^ ~ 2(a + s)(a" +6") ' 

inserting a = 1 gives the present day value, i.e.: 

^ {y - l)a + zs[ay- 2(1 + v)] + (/3 - 4)(1 + y) 
2(2 + z,)(l + y) 

with y — {\ + Zb)^" . Some straightforward considerations allows to derive the following constraints on q^r^: 



(10) 



(11) 
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(12) 
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It is convenient to solve Eq. (|lip with respect to in order to express this one as a function of qq and z^. It is : 



Zb 



a{l + Zs)+l3-{2 + z,){2qo + 2) 



-1 l/a 



a- P+{2 + z.s)i2qo + 2) 



- 1 



(13) 



Such a parametric density model is fully characterized by five parameters which are chosen to be the two asymptotic 
slopes {a, f3), the present day values of the deceleration parameter and of the Hubble constant (go, Hq) and the scaling 
redshift z^. As ini^, we will set (a,/?) = (3,4) and Zg = 3454 so that (goj-^^o) will be the two parameters to be 
constrained by the data. Any generalized EoS approach can be reduced to this scheme which is useful to fit the data. 
However, the phenomenological parameters, can often have a fundamental physics counterpart. 



C. Extended Gravity: the case of curvature quintessence 

There is no a priori reason to restrict the gravity action to a linear function of Ricci scalar R as in the General 
Relativity. Actually, higher order curvature invariants always emerge in quantum field theory formulated on curved 
spacetimes so it is worth considering the effect of such generalizations also on the late evolution of the universe. 

Replacing the Ricci scalar R with a (up to now) generic function /(i?) of Ricci scalar in the gravity lagrangian, the 
resulting field equations may still be recast in the Friedmann - like form provided that the total energy density and 
the total pressure are written as : 

Ptot = PM + Pcurv , Ptot — PM + Pcurv 

where the energy density and the pressure due to the higher - order curvature invariants are^^ : 



^ ^hfiR)-Rnm-iHRf"{R) 



f'{R) {2 



(14) 



{{2HR + R)riR) + R'riR) 



(15) 



This approach is particularly useful since it allows to interpret the non-Einstein part of gravitational interaction as a 
new fluid with energy density and pressure given by p4p and (|15p respectively. The barotropic factor for such a fluid 
turns out to be : 



-1 + 



nR)R+ 


'r"{R)R- 


- Hf"{R) 


IfiR) - Rf'iR)] 12 - 


3Hf"{R) 



(16) 



A leading role is played by the form chosen for f{R). As a simple choice^i, we set f{R) — foR" with /o constant. It 
is possible to find power law solutions of the Friedmann equations as a{t) = {t/to)" with: 



a 



2n^ - 3n + 1 
2-n 



(17) 



This gives a class of cosmological models with constant deceleration parameter : 



lit) = go = 



1 - a 2n^ - 2n - 1 



(18) 



2n2-3n+l ' 

In order to get an accelerated expansion (a > and go < 0) , the parameter n has to satisfy the following constraint : 
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1 - V3 



(19) 



I will refer to models with n in the first (second) range as CurvDown (CurvUp) models respectively and use the 
method described in the following section to constrain the two parameters (n, Hq) which completely assign the model. 
A comment is necessary at this point. Actually a simple power-law /(i?)-gravity model is not sufhcient to realistically 
reproduce a matter-dominated transient phase (needed for large scale structure formation) evolving into the today 
observed accelerated phase^i. However, also if more realistic models are needed'^^, it is useful to show that accelerated 
behaviour can be recovered simply by extending General Relativity without introducing further exotic ingredients 
(see also the contributions by Francaviglia, Nojiri, Odintsov and Schmidt in these Lecture Notes). 



III. METHODS TO CONSTRAIN MODELS 



In this section, I will discuss how cosmological models can be constrained, in principle, using suitable distance 
and/or time indicators. As a general remark, solutions coming from cosmological models have to be matched with 
observations by using the redshift z as the natural time variable for the Hubble parameter, i.e. 



H{z) 



1 



(20) 



Interesting ranges for z are: 100 < z < 1000 for early universe (CMBR data), 10 < z < 100 (LSS), < z < 10 (SNela, 
radio-galaxies, etc.). The method consists in building up a reasonable patchwork of data coming from different 
epochs and then matching them with the same cosmological solution ranging, in principle, from inflation to present 
accelerated era. In order to constrain the parameters characterizing the cosmological solution, a reasonable approach 
is to maximize the following likelihood function : 



C cx exp 



xHp) 



(21) 



where p are the parameters characterizing the cosmological solution. The merit function can be defined as : 



x^P) = E 

1=1 

Terms entering Ea. (j22p can be characterized as follows. For example, the dimensionless coordinate distances y to 
objects at redshifts z are considered in the first term. They are defined as : 



y {zi,p) - y, 



obs 



7^(p) - 1.716 
0.062 



^(p) - 0.469 
0.017 



(22) 



where E(z) = H{z)/Hq is the normalized Hubble parameter. This is the main quantity which allows to compare the 
theoretical results with data. The function y is related to the luminosity distance = (1 + z)r{z). A sample of data 
on y{z) for the 157 SNela is discussed in the Riess et al.- Gold dataset and 20 radio-galaxies are in'^^. These authors 
fit with good accuracy the linear Hubble law at low redshift {z < 0.1) obtaining the Hubble dimensionless parameter 
h — 0.664±0.008 . Such a number can be consistently taken into account at low redshift. This value is in agreement 
with Hq — 72±8 km s~^ Mpc~^ given by the HST Key projecli^^ based on the local distance ladder and estimates 
coming from time delays in multiply imaged quasars^ and Sunyaev- Zel'dovich effect in X-ray emitting clusters^. 
The second term in Eq. (f22|) allows to extend the z-range to probe y{z) up to the last scattering surface (z > 1000). 
The shift parameter^^^^^ TZ = \/Q.My{zis) can be determined from the CMBR anisotropy spectrum, where zis is the 
redshift of the last scattering surface which can be approximated as zis = 1048 (l -I- 0.00124a;j^'^'^^^) (1 + gi^%) with 
LUi = Vtih? (with i = b, M for baryons and total matter respectively) and (51,52) given in^^. The parameter ujh is 
constrained by the baryogenesis calculations contrasted to the observed abundances of primordial elements. Using 
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this method, the value uji, = 0.0214±0.0020 is found^. In any case, it is worth noting that the exact value of zis 
has a negligible impact on the results and setting z;^ = 1100 does not change constraints and priors on the other 
parameters of the given model. The third term in the function takes into account the acoustic peak of the large 
scale correlation function at 100 Mpc separation, detected by using 46748 luminous red galaxies (LRG) selected 
from the SDSS Main Sample^i^. The quantity 



is related to the position of acoustic peak where zlug = 0.35 is the effective redshift of the above sample. The 
parameter A depends on the dimensionless coordinate distance (and thus on the integrated expansion rate), on ^Im 
and E{z). This dependence removes some of the degeneracies intrinsic in distance fitting methods. Due to this reason, 
it is particularly interesting to include ^ as a further constraint on the model parameters using its measured value 
A = 0.469±0.017— . Note that, although similar to the usual introduced in statistics, the reduced (i-e-j the 
ratio between the and the number of degrees of freedom) is not forced to be 1 for the best fit model because of the 
presence of the priors on TZ and A and since the uncertainties ai are not Gaussian distributed, but take care of both 
statistical errors and systematic uncertainties. With the definition (|2ip of the likelihood function, the best fit model 
parameters are those that maximize C{p). 

Using the method sketched above, the classes of models studied here can be constrained and selected by observations. 
However, most of the tests recently used to constrain cosmological parameters (such as the SNela Hubble diagram and 
the angular size - redshift) are essentially distance - based methods. The proposal of Dalai et al.^^ to use the lookback 
time to high redshift objects is thus particularly interesting since it relies on a completely different observable. The 
lookback time is defined as the difi'erence between the present day age of the universe and its age at redshift z and 
may be computed as : 



where tn = ^/Hq = 9.78h~^ Gyr is the Hubble time (with h the Hubble constant in units of 100 km s~^ Mpc^^), and, 
as above, E{z,p) = H{z)/Hq is the dimensionless Hubble parameter and {p} the set of parameters characterizing 
the cosmological model. It is worth noting that, by definition, the lookback time is not sensible to the present day 
age of the universe so that it is (at least in principle) possible that a model fits well the data on the lookback time, 
but nonetheless it predicts a completely wrong value for to. This latter parameter can be evaluated from Eq. ((25|) by 
simply changing the upper integration limit from z to infinity. This shows that it is indeed a different quantity since 
it depends on the full evolution of the universe and not only on how the universe evolves from the redshift z to now. 
That is why this quantity can be explicitly introduced as a further constraint. As an example, let us now discuss how 
to use the lookback time and the age of the universe to test a given cosmological model. To this end, let us consider 
an object i at redshift z and denote by ti{z) its age defined as the difference between the age of the universe when 
the object was born, i.e. at the formation redshift zp, and the one at z. It is : 



where, in the last row, we have used the definition (|25p of the lookback time. Suppose now we have N objects and 
we have been able to estimate the age ti of the object at redshift for i = 1, 2, . . . , N. Using the previous relation, 
we can estimate the lookback time t°^^{zi) as : 




(24) 




(25) 




(26) 




[tT-u{z)]-K 

t"'-' - t,{z) - df , 



.obs 



= t\ 



(27) 



where f^" is the today estimated age of the universe and a delay factor can be defined as : 
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(28) 



The delay factor is introduced to take into account our ignorance of the formation redshift zp of the object. Actually, 
what can be measured is the age ti{z) of the object at redshift z. To estimate zp, one should use Eg. ([26)1 assuming 
a background cosmological model. Since our aim is to determine what is the background cosmological model, it is 
clear that we cannot infer zp from the measured age so that this quantity is completely undetermined. It is worth 
stressing that, in principle, df should be different for each object in the sample unless there is a theoretical reason to 
assume the same redshift at the formation of all the objects. If this is indeed the case (as I will assume later), then 
it is computationally convenient to consider df rather than zp as the unknown parameter to be determined from the 
data. Again a likelihood function can be defined as : 



At(p,d/)cxexp[-x?,(p,d/)/2] 



(29) 



with : 



Xu = 



1 



N-Np + l 



N 

E 



t*^'=°'-(z„p)-if^(z,) 



(30) 



where Np is the number of parameters of the model, at is the uncertainty on , ai the one on t°^''{zi) and the 
superscript theor denotes the predicted values of a given quantity. Note that the delay factor enters the definition 
of xft since it determines t'j^''{zi) from ti{z) in virtue of Eg. ip?)) . but the theoretical lookback time does not depend 
on df. In principle, such a method should work efficiently to discriminate among the various dark energy models. 
Actually, this is not exactly the case due to the paucity of the available data which leads to large uncertainties on the 
estimated parameters. In order to partially alleviate this problem, it is convenient to add further constraints on the 
models by using Gaussian priors^^ on the Hubble constant, i.e. redefining the likelihood function as : 



£(p) oc /;;t(p)exp 



h~h 



ohs 



cx exp[-x^(p)/2] 



(31) 



where we have absorbed df in the set of parameters p and have defined : 



X 



h-h 



obs 



(32) 



with h°^^ the estimated value of h and Uh its uncertainty. The HST Key project results^S can be used setting 
{h,afi) = (0.72,0.08). Note that this estimate is independent of the cosmological model since it has been obtained 
from local distance ladder methods. The best fit model parameters p may be obtained by maximizing £(p) which is 
equivalent to minimize the defined in Ea. (j32p . It is worth stressing that such a function should not be considered 
as a statistical in the sense that it is not forced to be of order 1 for the best fit model to consider a fit as successful. 
Actually, such an interpretation is not possible since the errors on the measured quantities (both ti and to) are not 
Gaussian distributed and, moreover, there are uncontrolled systematic uncertainties that may also dominate the error 
budget. Nonetheless, a qualitative comparison among different models may be obtained by comparing the values of 
this pseudo even if this should not be considered as a definitive evidence against a given model. Having more 
than one parameter, one obtains the best fit value of each single parameter pi as the value which maximizes the 
marginalized likelihood for that parameter defined as : 



£p. cx 



dpi 



dpn C{p) 



(33) 



After having normalized the marginalized likelihood to 1 at maximum, one computes the la and 2a confidence limits 
(CL) on that parameter by solving Cp. = exp (—0.5) and Cp. = exp (—2) respectively. In summary, taking into 
account the above procedures for distance and time measurements, one can reasonably constrain a given cosmological 
model. In any case, the main and obvious issue is to have at disposal sufficient and good quality data sets. 
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Color age 


Scatter age 


z 


N 


Age (Gyr) 


Ref 


z 


N 


Age (Gyr) 


Ref 


0.60 


1 


4.53 


32 


0.10 


55 


10.65 


34 


0.70 


3 


3.93 


32 


0.25 


103 


8.89 


34 


0.80 


2 


3.41 


32 


1.27 


1 


1.60 


35 



TABLE I: Main properties of the cluster sample used for the analysis. The data in the left part of the Table refers to clusters 
whose age has been estimated from the color of the reddest galaxies (color age), while that of clusters in the right part has 
been obtained by the color scatter (scatter age). For each data point, I give the redshift z, the number A'^ of clusters used, the 
age estimate and the relevant reference. 

IV. SAMPLES OF DATA TO CONSTRAIN MODELS: THE CASE OF LSS 

In order to apply the method outlined above, we need a set of distant objects whose age can be somehow estimated. 
For example, clusters of galaxies seem to be ideal candidates in this sense since they can be detected up to high redshift 
and their redshift, at formation epoch^" is almost the same for all the clusters. Furthermore, it is relatively easy to 
estimate their age from photometric data only. To this end, the color of their component galaxies, in particular the 
reddest ones, is needed. Actually, the stellar populations of the reddest galaxies become redder and redder as they 
evolve. It is just a matter, then, to assume a stellar population synthesis model, and to look at how old the latest 
episode of star formation should be happened in the galaxy past to produce colors as red as the observed ones. This 
is what is referred to as color age. The main limitation of the method relies in the stellar population synthesis model, 
and on a few (unknown) ingredients (among which the metallicity and the star formation rate). The choice of the 
evolutionary model is a key step in the estimate of the color age and the main source of uncertainty"^*^. An alternative 
and more robust route to cluster age is to consider the color scatter (see^ for an early application of this approach). 
The argument, qualitatively, goes in this way: if galaxies have an extreme similarity in their color and nothing is 
conspiring to make the color scatter surreptitiously small, then the latest episode of star formation should happen in 
the galaxy far past, otherwise the observed color scatter would be larger. Quantitatively, the scatter in color should 
thus be equal to the derivative of color with time multiplied the scatter of star formation times. The first quantity 
may be predicted using population synthesis models and turns out to be almost the same for all the evolutionary 
models thus significantly reducing the systematic uncertainty. We will refer to the age estimated by this method as 
scatter age. The dataset we need to apply the method described in the previous section may now be obtained using 
the following procedure. First, for a given redshift Zi, we collect the colors of the redddest galaxies in a cluster at that 
redshift and then use one of the two methods outlined above to determine the color or the scatter age of the cluster. 
If more than one cluster is available at that redshift, we average the results from different clusters in order to reduce 
systematic errors. Having thus obtained ti{zi), we then use Ea.(P7|) to estimate the value of the lookback time at that 
redshift. Actually, what we measure is t°^''{zi) + df that is the quantity that enters the definition l[30|) of xit ^'^^ then 
the likelihood function. To estimate the color age, following^, it is possible to choose, among the various available 
stellar population synthesis models, the Kodama and Arimoto one^, which, unlike other models, allows a chemical 
evolution neglected elsewhere. This gives us three points on the diagram z vs. t°j^^ obtained by applying the method 
to a set of six clusters at three different redshifts as detailed in Table 1. Using a large sample of low redshift SDSS 
clusters, it is possible to evaluate the scatter age for clusters age at 2 = 0.10 and z — 0.25'^'^. Blakeslee et ah-- applied 
the same method to a single, high redshift {z — 1.27) cluster. Collecting the data using both the color age and the 
scatter age, we end up with a sample of ^ 160 clusters at six redshifts (listed in Table 1) which probe the redshift 
range (0.10, 1.27). This nicely overlaps the one probed by SNela Hubble diagram so that a comparison among our 
results and those from SNela is possible. I assume a cr = 1 Gyr as uncertainty on the cluster age, no matter what is 
the method used to get that estimate. Note that this is a very conservative choice. Actually, if the error on the age 
were so large, the color - magitude relation for reddest cluster galaxies should have a large scatter that is not observed. 
I have, however, chosen such a large error to take qualitatively into account the systematic uncertainties related to 
the choice of the evolutionary model. 

Finally, we need an estimate of tg''^ to apply the method. Following Rebolo et ali^, one can choose {t'^'' ,at) = 
(14.4, 1.4) Gyr as obtained by a combined analysis of the WMAP and VSA data on the CMBR anisotropy spectrum 
and SDSS galaxy clustering. Actually, this estimate is model dependent since Rebolo et aL?^ implicitly assumes that 
the ACDM model is the correct one. However, this value is in perfect agreement with = 12.6^2 4 Gyr determined 
from globular clusters age-^ and > 12.5 ± 3.5 Gyr from radioisotopes studies'"^. For this reason, one is confident 
that no systematic error is induced on the adopted method using the Rebolo et al. estimate for ig'"' even when testing 
cosmological models other than the ACDM one. 
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FIG. 2: Comparison among predicted and observed values of t = tL{z) + df for the best fit QCDM model. 

V. TESTING THE DIFFERENT CLASSES OF COSMOLOGICAL MODELS 

The method outlined above can be applied to the dark energy models described in Sect. 2 in order to constrain their 
parameters and see if they are viable candidates to explain cosmic acceleration. To this aim, let us note that each one 
of the presented models is roughly described by few parameters which are^ {^M,h, w) for the QCDM model, (qoj h) 
for the parametric density model and (n, h) for the curvature quintessence. For all the three classes of models, there 
is still another parameter entering the fit, that is the delay factor d/, which we will marginalize over since it is not 
interesting for our aims. Let us first consider the QCDM model. The main results are plotted in Figs. [2] and [3l In the 
first plot, we compare the estimated clusters age with the quantity: 

T{z)=tL{z) + df (34) 
using the best fit values for the model parameters and the delay factor which turn out to be : 

= (0.25,0.70,-0.81) , d/ = 4.5Gyr 

giving ~ 0.04. The value for the best fit parameters (both for the QCDM model and the other dark energy 
models considered) turns out to be quite small suggesting that errors have been seriously overestimated. This is not 
surprising given the arbitrary way we have fixed the uncertainties on the estimated age of the clusters. That this is 
likely to be the case is also suggested by a qualitative argument. One could rescale the errors on tf'{zi) in such a way 
that \^ — 1 for the best fit model. Since for the best fit QCDM model, ~ x^f ^^^^ leads to multiply by almost 1/5 
the uncertainties on t°^''{zi). If the error on to were negligible, this means that we should reduce the uncertainty on 
the cluster age from 1 Gyr to 0.2 Gyr that is indeed a more realistic value. The presence of an error on t'^'^ slightly 
complicates this argument, but does not change the main conclusion. We are thus confident that the very low value of 
the g6t for the best fit model is only due to overestimating the uncertainties on the clusters ages. However, we 

do not perform any rescaling of the uncertainties since, to this end, we should select a priori a model as the confidence 
one which is contrary to the adopted philosophy by which we want to point out the degeneracy. It is worth stressing, 
however, that such rescaling does not affect anyway the main results. 

Fig. [3] shows the Icr and 2<t confidence levels in the {^Miw) plane obtained by marginalizing over the Hubble 
constant and the delay factor. Two interesting considerations may be drawn from that plot. First, we note that also 
phantom models (i.e. models with w < — 1 violating the weak energy condition) are allowed by the data. This is in 
agreement with recent results coming from fitting the QCDM model to the SNela Hubble diagram and the CMBR 
anisotropy spectrum^. Unfortunately, a direct comparison is not possible since the marginalized likelihood is too flat 
to get any constraints on w so that all the values in the range tested (— 2 < w < 1/3) are well within the la CL. 
Although essentially due to the paucity of the data, this result is also a consequence of having not used any prior 
on U,M as it is usually done in other analyses. By using the procedure described at the end of Sect. 3, we obtain the 
following estimates for the other QCDM parameters : 

VLm e (0.13,0.39) , h e (0.63,0.77) (la CL) , 




e (0.01,0.62) , /i e (0.56,0.84) (2cr CL) 
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FIG. 3: The la and 2a confidence regions in the plane {Qm,w) for the QCDM model after marginalizing over the Hubble 
constant h and the delay factor df. 
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FIG. 4: Comparison among predicted and observed values of r = tL{z) + df for the best fit ACDM model. 

Given that we are able to give only weak constraints on the QCDM model, from now on I will dedicate attention to 
the case w — ~1, i.e. to the ACDM model and do not discuss anymore the results for the QCDM model. The best 
fit parameters for the cosmological constant model turn out to be : 

(r^M, h) = (0.22, 0.71) , df = 4.05 Gyr (x^ ~ 0.09) 

that gives rise to the curve r(z) shown in Fig. 31 while Fig. O reports the confidence regions in the {^Im, h) plane after 
marginalizing over the delay factor. From the marginalized likelihood functions, we get : 

VLm e (0.10,0.35) , h e (0.63,0.78) {la CL) , 



e (0.06,0.59) , h e (0.56,0.85) {2a CL) . 

The ACDM model has been widely tested against a large set of different astrophysical data. This offers us the 
possibility to cross check both the model and the validity of method. It is instructive, in this sense, to compare our 
results with those coming from the fit to the SNela Hubble diagram. For instance, Barris et al.^ used a set of 120 
SNela up to z = 1.03 finding Qm — 0.33 as best fit value with a large uncertainty (not quoted explicitly, but easy to 
see in their Fig. 12) in good agreement with our result. A more recent result has been obtained by Riess et alii which 
have used a SNela Hubble diagram extending up to z = 1.55 and have found = 0-291qq3 still in agreement with 
our result. The ACDM model has also been tested by means of the angular size - redshift relation. Using a catalog of 
ultracompact radio sources and taking carefully into account systematic uncertainties and selection effects, Jackson^ 
has found Qm = 0.24tg °? in almost perfect agreement with our estimate. 
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FIG. 5: The la and 2a confidence regions in the plane (rinijh) for the ACDM model. 
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FIG. 6: Comparison among predicted and observed values of t — tL{z) + df for the best fit parametric density model. 

Neither Barris et al."^ nor Riess et alX quote a best fit value for h since this parameter is infinitely degenerate with 
the supernovae absolute magnitude M when fitting to the SNela Hubble diagram. Nonetheless, SNela may still be 
used to determine h by fitting the linear Hubble law to the low redshidt (z < 0.1) SNela. Using this method, Daly & 
Djorgovski'^^ have found h = 0.664 ± 0.08 in good agreement with the presented result. Moreover, it is worth noting 
that our estimate for h turns out to be in agreement with estimates coming from different methods such as various local 
standard candles^^, the Sunyaev- Zel'dovich effect in galaxy clusters^ and time delays in multiply imaged quasars^. 
Finally, let us quote the results obtained by Tegmark et aL— which have performed a combined fit of the ACDM model 
to both the WMAP data on the CMBR anisotropy spectrum and the galaxy power spectrum measured by more than 
200,000 galaxies surveyed by the SDSS collaboration. They find Qm = 0.30 ± 0.04 and h = O.TOto ^ot in very good 
agreement with our results. Perhaps, the most interesting result of testing the ACDM model with lookback method is 
not the success of this model (since it has yet been shown by a lot of evidences) , but the substantial agreement among 
these estimates of the parameters (fi^/, h) and those coming from distance measurements. This is quite encouraging 
since it is an important successful cross check which, in general, could be applied to any cosmological models. 

Let us now examine the results obtained for the parametric density model, pictorially shown in Figs. [6] and [T] The 
best fit model is obtained for the following values of the parameters {qq, h) and of the delay factor : 

(go, h) = (-0.68, 0.71) , df = 4.20 Gyr (x' ^ 0.07) . 
Marginalizing over df, we get : 

qo e (-0.81,-0.47) , h e (0.64,0.78) (Icr CL) , 




go e (-0.89,-0.24) , /i e (0.58,0.85) (2cr CL) 
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FIG. 7: The la and 2a confidence regions in the plane (go, h) for the parametric density model. 
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FIG. 8: The la and 2a confidence regions in the plane (n, h) for curvature quintessence in the CurvUp regime. 

Note that the 2a CL on the go parameter has been truncated at the upper end since it formally extends to values 
higher than the physically acceptable one. 

In^'^, the parameters of this model have been constrained by using both the SNela Hubble diagram and the angular 
size - redshift relation. In particular, fitting the model to the SNela Hubble diagram gives h = 0.64lQ Qg, while the 
physically acceptable range for go turns out to be in agreement with the data for go = —0.42 as best fit value. The 
present day deceleration parameter go is better constrained using the data listed in Jackson^- to perform the angular 
size - redshift test thus obtaining go — — 0.641q 'ijj^^. Both these results are in very good agreement with the present 
estimates so that we conclude that the parametric density model is a viable candidate alternative to the ACDM model 
but this is a signal of degeneracy among the two models. 

Finally, let us discuss the results for the curvature quintessence. Figs.[S] and [S] report the confidence regions in the 
plane (n, h) for the CurvUp and CurvDown regime respectively after marginalizing over the delay factor. The first 
striking feature is that the contour plots are not closed so that the marginalized likelihood function gives only an 
upper (lower) limit to the parameter n in the CurvUp (CurvDown) regime. Formally, we get the following estimates 
for the best fit values in the CurvUp and CurvDown regime respectively : 

(n, h) = (1.367, 0.71) , df = 4.80 Gyr (x^ 0.23) , 





(n, h) = (-0.367, 0.74) , df = 4.80 Gyr (x^ ^ 0.21) , 
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FIG. 9: The la and 2a confidence regions in the plane (n, h) for curvature quintessence in the CurvDown regime. 

but the best fit value for n actually lies outside the investigated range for this parameter. Being the confidence regions 
open, it is meaningless to give constraints on h, but nonetheless it is possible to infer the following limits : 

n < 1.402 at la CL , n< 1.424 at 2a CL , 



n > -0.508 at Icr CL , n> -0.606 at 2a CL , 

where the first row refers to the CurvUp regime and the second to the CurvDown one. These limits do not contradict 
the ranges determined by fitting to the SNela Hubble diagrami^, but seems to be quite non-realistic. Actually, the 
deceleration parameter corresponding to these values of n is quite small (qq ~ —0.01) thus contradicting the evidence 
in favor of an accelerating universe. Moreover, the results in^^ have been obtained by using an old sample of SNela, 
including some SNela that have now been discarded from the Gold set in Riess et alJ. On the other hand, it is 
worth noting that fitting a power law scale factor to the angular size - redshift relation for compact radio sources gives 
a ~ 1— which, by using Eg. (fT7|) . translates in an estimate for n that is in good agreement with the presented result. 

VI. DISCUSSION AND CONCLUSIONS 

The impressive amount of data indicating a spatially flat universe in accelerated expansion has posed the problem 
of dark energy and stimulated the search for cosmological models able to explain such unexpected behavior. Many 
rival theories have been proposed to solve the puzzle of the nature of dark energy ranging from a rolling scalar field 
to a unified picture where a single exotic fluid accounts for the whole dark sector (dark matter and dark energy). 
Moreover, modifications of the gravity action has also been advocated. Although deeply different in their underlying 
physics, all these scenarios share the common feature of well reproducing the available astrophysical data giving 
rise to a degeneracy problem which is easily showed also in the straightforward analysis presented here. It is worth 
stressing, however, that the most widely used cosmological tests and in particular the SNela Hubble diagram and 
the angular size - redhisft relation) are essentially based on distance measurements to high redshift objects and are 
thus affected by similar systematic errors. It is hence particularly interesting to look for methods which are related 
to the estimates of different quantities. Being affected by other kinds of observational problems, such methods could 
be considered as cross checks for the results obtained by the usual tests and they should represent complementary 
probes for cosmological models. The technique we have devised here is a first step in this direction. We have used the 
present age of the universe and the lookback time to galaxy clusters to build a sort of time diagram of the universe 
in order to reconstruct its age evolution. Relying on stellar evolutionary codes, the estimate of the lookback time is 
related to a different astrophysics than the distance based methods and it is thus free from any problem connected 
with the evolution of standard candles (such as the SNela absolute magnitude and the intrinsic linear size of radio 
sources). Actually, this technique could be affected by its own systematics (such as, for instance, the degeneracy 
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between age and metallicity), but these may be better controlled. Moreover, comparing the results thus obtained 
with those coming from distance based methods allows to strengthen the conclusions suggested by both techniques. 
Motivated by these considerations, I have first applied the method to the ACDM model in order to estimate the 
present day values of the matter density parameter Qm and of the dimensionlcss Hubble constant h obtaining (at the 
la CL) : 



CIm — 0.221q :[q , h — 0.7lto;gg . 

These values are in agreement with the previous estimates in literature. This can be considered an independent 
confirmation not only of the viability of the ACDM model, but also of the method. It is worth noting that the ACDM 
scenario receives further support from this test and henceforth the cosmological constant A still remains the best 
candidate to explain the dark energy puzzle from an observational point of view. Nonetheless, the ACDM model is 
severely affected by conceptual problems so that it is worthwhile to look at alternative approaches. This has stimulated 
a plenty of models where the cosmic acceleration is due to a dominant scalar field rolling down its potential. However, 
such a scheme, on the one hand, still does not solve the coincidence problem and, on the other hand, is plagued by 
the unidentified nature of the scalar field itself and the ignorance of its self interaction potential. These considerations 
have opened the way to different models that are able to give an accelerated expansion without the need of scalar 
fields. Two of these approaches have been tested in these lectures. 

Ir>i^, phenomenological unified models have been proposed where a single fluid with a given energy density assigned 
by few parameters is able to fit both the SNela Hubble diagram and the angular size - redshift relation for ultracompact 
radio sources. These parametric density models have been tested here by lookback time method and the following 
estimates (at the la CL) for its characterizing parameters have been obtained: 



qo^~0.68t^Qil , = 0.71 ± 0.07 , 

in agreement with previous testsi^,. It is worth stressing that such a method turns out to be more efficient than the 
usual fit to the SNela Hubble diagram for the parametric density models. It is also interesting to note that both go 
and to predicted by these models are in almost perfect agreement with those computed for the ACDM model with the 
parameters (D,M,h) as discussed above. This is not surprising since Eq.® shows that, nowadays, the energy density 
of the parametric density model is very similar to that of the ACDM model so that these two scenarios share most 
of the observable properties referred to today quantities. However, this does by no way mean that the two models 
are the same. Actually, it is the underlying philosophy that is completely different having now a single fluid rather 
than a cosmological constant dominating over the matter term. Both models predict similar values of the today 
observed quantities simply because they are tied to reproduce the same data and not because they are two different 
manifestations of the same underlying physics. Another possible approach to the cosmic acceleration is to consider 
dark energy and dark matter as signals of the breakdown of the Einstein General Relativity at some characteristic 
scalc^^. In this picture, the universe is still dominated by standard matter, but the Friedmann equations have to be 
modified as consequence of a different gravity action. This philosophy inspired curvature quintessence scenarios where 
an effective dark energy is related to the properties of the function f{R) which replaces the Ricci scalar R in the 
gravity Lagrangianl^. The coupling between matter and curvature for a theory with f{R) = foR" gives rise to fourth 
order nonlinear differential equations for the evolution of the scale factor. Power law solutions, a ~ i", are possible, 
provided that a is linked to n by Eq. (|17[) . Such models are particularly attractive from a theoretical viewpoint (since 
they allow to give a purely geometric interpretation of the dark energy) and has also been shown to fit well SNela 
Hubble diagram^. Here it is shown that, while they are able to pass successfully the lookback time tests, only weak 
constraints can be imposed on the parameter n with the best fit value lying in a region corresponding to decelerating 
rather than accelerating models. We are thus tempted to conclude that this scenario could be rejected with some 
confidence. Actually, this conclusion is also suggested by the results in Riess et ali^ which, using the SNela Hubble 
diagram up to z = 1.55, have detected a change in the sign of the deceleration parameter q. The power law solutions 
predict, on the contrary, a constant q so that they are ruled out by the result in Riess et al.^. It is important to stress, 
however, that this is not a definitive exclusion of the curvature quintessence scenario. Actually, I have only considered 
simple power law solutions. It is indeed conceivable that considering more realistic models'^'', the evolution of the 
scale factor gives a past deceleration followed by the present acceleration thus leading the curvature quintessence 
models in agreement with what is suggested by the extended SNela Hubble diagram and other tests^. Having tested 
three different models, it is worth asking what is the better one. Unfortunately, this is not possible on the basis of 
the test results only. Comparing the values, one could naively conclude that the parametric density model is the 
best one since it gives the lowest x^. However, while the difference among the values is significant between the 
parametric density model and the curvature quintessence scenario (almost an order of magnitude), it is too small 
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(0.07 vs. 0.09) to conclude that the ACDM model is disfavored. We have thus to conclude that this test alone is not 
able to discriminate among these two dark energy candidates. On the other hand, Sandvik et al4^ recently claimed 
that UDE models are not viable because the growth of density perturbations will lead to matter power spectrum in 
disagreement with what is observed. This should be an evidence against the parametric density model. However, it is 
worth noting that Sandvik et al. explicitely consider the generalized Chaplygin gas model which is characterized by 
a negative squared sound speed that seems to be the main cause of the anomalous growth of perturbations. For the 
parametric density model, the sound speed is always positive definite so that it is likely that the argument of Sandvik 
et al. should be at least reconsidered. Finally, I would briefly comment on the possibility to ameliorate methods by 
increasing the maximum redshift probed. To this aim, galaxy clusters do not appear as good candidates since it is 
quite difficult to detect a significative number of member galaxies up to redshifts larger than ~ 1.3. However, high 
redshift galaxies may be taken into account provided that they are detected in as more photometric bands as possible. 
This latter requirement is fundamental since it allows not only to better estimate the age of the galaxy, but also to 
infer constraints on its formation redshift zp that cannot be assumed to be the same for all the galaxies (as it has 
been possible for clusters). Some techniques allow to discover galaxies up to redshift z ~ 6^ and thus let us hope to 
measure lookback time up to such high redshift. With its ability of both furnishing multicolor photometry of high 
redshift galaxies (and thus better estimates of their color ages and formation redshift) and increasing the number of 
SNela, surveys like the GOODS one'*- seems to be the most promising source of cosmological constraints in the near 
future. In any case, from one hand we need some experimentum crucis capable of removing the degeneracy for a 
reasonably large redshift range and, from the theoretical viewpoint, we need a physically reliable cosmological model, 
emerging from some fundamental theory, without the conceptual shortcomings of ACDM. 
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EoS for Equation of State. 

5* It is worth stressing that this model may be interpreted not only as comprising a single fluid with an exotic equation of 
state, but also as made of dark matter and scalar fleld dark energy or in the framework of modified Friedmann equations. 
For pedagogical reasons, I prefer the UDE interpretation even if the results do not depend on this assumption. 

5® The need for priors to reduce the parameter uncertainties is often advocated for cosmological tests. For instance, in Ref.— 
a strong prior on is introduced to constrain the dark energy equation of state. It is likely, that extending the dataset to 
higher redshifts and reducing the uncertainties on the age estimate will allow to avoid resorting to priors. 
It is worth stressing that, in literature, the cluster formation redshift is defined as the redshift at which the last episode of 
star formation happened. In this sense, we should modify our definition of df by adding a constant term which takes care of 
how long is the star formation process and what is the time elapsed from the beginning of the universe to the birth of the 
first cluster of galaxies. For this reason, it is still possible to consider the delay factor to be the same for all clusters, but 
it is not possible to infer zf from the fitted value of df because we do not know the detail of star formation history. This 
approach is particular useful since it allows to overcome the problem to consider lower limits of the universe age at z rather 
than the actual values. 

®^ We drop the subscript from SIa/.o since it does not give rise to any confusion here. Also, we use the dimensionless parameter 
h instead of the Hubble constant Hq. 



